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Abstract. An explicit determination of all local conservation laws of kinematic type on 
moving domains and moving surfaces is presented for the Euler equations of inviscid com¬ 
pressible fluid flow in curved Riemannian manifolds in n > 1 dimensions. All corresponding 
kinematic constants of motion are also determined, along with all Hamiltonian kinematic 
symmetries and kinematic Casimirs which arise from the Hamiltonian structure of the in¬ 
viscid compressible fluid equations. 


1. Introduction 

The study of topological, geometrical, and group-theoretic aspects of fluid equations in 
dimensions n > 1 has attracted considerable interest [H [5l [6] in the mathematical theory of 
fluid flow. Two central topics in studying these aspects are Hamiltonian structures EE] 
and conserved integrals PED]. 

For the Euler equations of inviscid compressible fluid flow in multi-dimensional flat man¬ 
ifolds M”, there is a fairly complete picture of conserved integrals that arise from local con¬ 
servation laws (i.e. continuity equations) of kinematic type and vorticity type on domains 
moving with the fluid. A kinematic conservation law, like mass, energy, momentum and an¬ 
gular momentum, refers to a continuity equation in which the conserved density and spatial 
flux involve only the fluid velocity, density and pressure, in addition to the time and space 
coordinates. In contrast, a vorticity conservation law, such as helicity in three dimensions as 
well as circulation and enstrophy in two dimensions, refers to a continuity equation where 
the conserved density and spatial flux have an essential dependence on the curl of the fluid 
velocity. These two classes of continuity equations comprise all of the local conservation laws 
found to-date for inviscid compressible fluid flow in M” (with n > 1). 

An explicit classihcation of kinematic and vorticity conservation laws on moving domains 
is known |2] in the case of inviscid compressible fluid flow with a barotropic equation of state 
for the pressure. In particular, the vorticity conservation laws are comprised by helicity 
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in all odd dimensions n > 3 and enstrophy in all even dimensions n > 2, while the only 
kinematic conservation laws apart from mass, energy, momentum and angular momentum 
consist of Galilean momentum which holds for all equations of state, plus a similarity energy 
and a dilational energy which arise for polytropic equations of state where the pressure is 
proportional to a special dimension-dependent power 7 = l-|-2/nof the density. Here a 
moving domain refers to a closed volume in M” that is transported along the streamlines of 
the fluid. 

A similar classihcation has been obtained recently [3] for inviscid non-isentropic compress¬ 
ible fluid flow in M” (with n > 1 ), where the entropy is conserved only along streamlines and 
the pressure is given by an equation of state in terms of both the fluid density and entropy. 
In this case, helicity and enstrophy are no longer conserved, but in all even dimensions n > 2 
there is a vorticity conservation law given by an entropy circulation (which vanishes when¬ 
ever the fluid is irrotational or isentropic), plus there is one extra kinematic conservation law 
consisting of volumetric entropy in any dimension. Both of these conservation laws hold for 
all equations of state. 

Much less is known, however, about conserved integrals for inviscid compressible fluid 
flow in multi-dimensional curved manifolds. One general result is that all of the vorticity 
conservation laws on moving domains for fluid flow in M"" have a natural generalization to 
curved Riemannian manifolds, because these conservation laws arise as Casimir invariants 
from the Hamiltonian structure of the Eulerian fluid equations [T 6 l [H [TT] . A related result [1] 
is that recently these conservation laws have been further generalized to lower-dimensional 
surfaces that move with the fluid, providing new conserved integrals on moving surfaces of 
any dimension in flat and curved manifolds. 

The present paper will settle the open question of explicitly determining all local conserva¬ 
tion laws of kinematic type on moving domains and moving surfaces for inviscid compressible 
fluid flow in curved Riemannian manifolds. In particular, any such conservation laws will 
be found that hold only for ( 1 ) special dimensions of the manifold or the surface; ( 2 ) special 
conditions on the geometry of the manifold or the surface; (3) special equations of state. 
Importantly, the general form of these kinematic conservation laws will be allowed to de¬ 
pend on the intrinsic Riemannian metric, volume form, and curvature tensor of the manifold 
or the surface. All kinematic constants of motion that arise from the resulting kinematic 
conservation laws also will be determined. 

A sequel paper will address the remaining open problem of determining whether the known 
local conservation laws of vorticity type on moving domains and moving surfaces are complete 
for inviscid compressible fluid flow in flat and curved manifolds. 

Note, in all of this work, the fluid is assumed to £11 the entire manifold. For results on 
conservation laws of fluid flow with a free boundary, see Ref. |14j. 

In section |2l first a summary of the Euler equations of inviscid compressible fluid flow 
in n-dimensional manifolds is given. Next the formulation of local conservation laws, con¬ 
served integrals, and constants of motion is discussed for general hydrodynamic systems in 
n-dimensional manifolds, and this formulation is adapted to moving domains and moving 
surfaces. Finally, necessary and sufficient determining equations are presented for directly 
finding all conserved densities of kinematic type on moving domains and moving surfaces for 
the Eulerian fluid equations. 
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The main results giving an explicit classification of all kinematic conserved densities on 
moving domains and moving surfaces for inviscid compressible fluid flow in u-dimensional 
manifolds are presented in section [3l A corresponding classification of kinematic constants of 
motion is also stated, along with Hamiltonian kinematic symmetries and kinematic Casimirs 
which arise from the Hamiltonian structure of the inviscid compressible fluid equations. 

The proof of these results is carried out in section m by solving the determining equations 
from section [21 The steps are carried out using tensorial index notation which is summarized 
in an Appendix. 

One interesting feature of the classifications is that special equations of state in which the 
pressure depends only on the entropy of the fluid are considered. For any such equation 
of state, new local conservation laws describing a generalized momentum and energy which 
depend on the entropy are found to arise for non-isentropic compressible fluid flow. 

Some concluding remarks are made in section O 

2. Preliminaries 

The Eulerian fluid equations in M"" are, in general, given 
mass density p, the entropy S, and the pressure P by 

Ut + u- Vu = —p“^VP, 

Pi + V • (pu) = 0, 

Fi + M-= 0, 

together with a general equation of state P = P(p, S). 

To generalize the Eulerian fluid equations to an n-dimensional manifold M [6], the only 
structure needed on M is a Riemannian metric g. Let V be the metric-compatible co¬ 
variant derivative determined by Vp = 0, and write grad and div for the contravariant 
gradient operator and the covariant divergence operator defined by ^JV = g{^,grad) and 
g{grad,^) = div ^ holding for an arbitrary vector field ^ on M. These operators are the 
natural Riemannian counterparts of the gradient V and divergence V- operators in R"'. For 
later use, let e be the volume form normalized with respect to p, and let e be the dual vol¬ 
ume tensor, satisfying Ve = 0 and p(e,e) = n\. Let Riem = [V, V] be the curvature tensor 
determined from p, and let R be the scalar curvature. Also, let Grad and Div denote the 
total contravariant gradient and the total covariant divergence, and let Dt denote the total 
time derivative, which are respectively defined by grad , div , dt acting via the chain rule. 

In this geometric notation, the covariant generalization of the fluid velocity equation fl2.1l) 
from R” to M is given by 

Ut +{u\V)u =—p~^ gradP (2.4) 

where u is the fluid velocity vector on M. Similarly the covariant equations for the fluid 
mass density p and entropy S on M are given by 

pt + div (pu) = 0, (2.5) 

St + u\VS = 0 . (2.6) 

A general equation of state is given by 

P = P(p, S) 
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in terms of the velocity u, the 

( 2 . 1 ) 

( 2 , 2 ) 

(2.3) 


(2.7) 




which closes the system fl2.4p - fl2.6p . Through a standard thermodynamic relation, the pres¬ 
sure P determines an associated internal (thermodynamic) energy which is defined by p!3lfT8] 

e(p,S) = j f,-^P(p,S)dp. (2.8) 

A transcription between geometric notation and tensorial index notation is provided in 
the beginning of Appendix B. 


2.1. Local conservation laws on moving domains. For any hydrodynamic system in 
a Riemannian manifold M, local conservation laws are described by a covariant continuity 
equation 

DtT + DivX = 0 (2.9) 

holding for all formal solutions of the system, where T and X are some functions of the hy¬ 
drodynamic variables and their spatial derivatives, as well as the time and space coordinates 
t, X. Physically, the scalar function T is a conserved density while the vector function X is a 
spatial flux. Note that, through their dependence on x, both T and X are allowed to depend 
on the metric tensor g, volume tensor e, and curvature tensor Riem. 

Consider any domain (i.e. an orientable closed spatial volume) V in M through which 
the fluid is flowing, and let z> be the outward unit normal vector on the boundary dV. In 
integral form on V, the continuity equation fl2.9p is equivalently given by 

^ [TdV = -[ g{X,u)dA (2.10) 

Jv Jav 

where dV = e is the volume n-form (dual of the volume tensor e), and dA = Dje is the 
hypersurface area n — 1-form in terms of the normal vector z/. 

A physically more useful form for expressing hydrodynamic conservation laws fl2.9p and 
fl2.10p is obtained by considering a domain V(t) that moves with the fluid. In particular, let 
each point x® G V{t) be transported along streamlines in the fluid, as defined by dx^/dt = 
CuP = nj Vx® {i = 1,..., n), where u is the fluid velocity vector and x® are local coordinates 
on M. Introduce the material (advective) derivative 

^t = Dt + Cu ( 2 . 11 ) 

where Cu denotes the Lie derivative with respect to the fluid velocity u. Then the spatial 
flux through the moving boundary dV{t) is given by 

$ = X-Tm (2.12) 

which is related to the conserved density T by the transport equation 

T)tT =—{div u)T — Div ^ (2-13) 


where 


div u = —e\Cu^ 
n\ 


(2.14) 


represents the expansion or contraction of an infinitesimal volume moving with the fluid. 
The corresponding integral form of the transport equation (I2.13p in V{t) is expressed as [1] 

d 


-r / TdV = 


dt 


'V{t) 


iav{t) 




(2.15) 








which is called a conserved integral on a moving domain in the fluid. As shown by equation 
fl2.15p . the integral expression TdV will be a constant of motion on V(t) C M if the net 
flux across the domain boundary dV(t) vanishes. 

Both the conserved integral fl2.15p and the underlying transport equation fl2.13p have an 
alternative formulation using differential forms, which generalizes in a simple way to moving 
surfaces. The following transport identity will be needed. Let S(t) C M be an orientable 
p-dimensional submanifold transported along the fluid streamlines, with 1 < p < n. Then 
for any p-form a, 

f ct = f Dtcy. (2-16) 

ls{t) Js{t) 


dt 


holds identically. (See Ref.[T] for a proof). 

This identity can be applied to the volume integral in equation fl2.15p . while the hypersur¬ 
face integral in equation fl2.15p can be converted into a volume integral by Stokes’ theorem, 
yielding 


dt 


[ TdV+ [ g{^,u)dA= [ Di(Te)+ d(<f)Je). 

Jv{t) JdV{t) Jv(t) 


(2.17) 


This integral equation holds on an arbitrary moving domain V(f) iff the integrand n-form 
vanishes. Hence the density T and flux $ satisfy 


2)t(Te) + d(<f)Je) = 0 


(2.18) 


which is equivalent to the transport equation (12.131) expressed in terms of differential forms. 
Note that here d is the exterior derivative acting as a total (spatial) derivative. 


2.2. Local conservation laws on moving surfaces. Let 1 < p < n — 1 and consider any 
p-dimensional surface (i.e. an orientable submanifold) S{t) in M that moves with the fluid, 
whereby each point x* G S(t) is transported along the fluid streamlines, dx^/dt = LuX^ = 
mJ Vx* (z = 1,..., n), in local coordinates on M. 

A conserved integral on a moving surface S{t) is the integral continuity equation [T] 


d 

dt 


OL = — 


'S{t) 


ias{t) 


/3 


(2.19) 


for a p-form density a and a p — 1-form flux (3 that are some functions of the hydrodynamic 
variables and their spatial derivatives, and the time and space coordinates t, x, holding for 
all formal solutions of the hydrodynamic system. The dependence of cx and /3 on x allows 
them to depend on any geometrical tensors defined on the surface S{t). 

The integral expression a will be a constant of motion on S(t) C M when the flux 
integral is zero for every formal solution of the hydrodynamic system. If S(t) is boundaryless 
then every conserved integral fl2.19p yields a constant of motion. Alternatively, if S(t) has a 
boundary dS{t) ^ 0 then a conserved integral fl2.19p yields a constant of motion only when 
(3 = d 7 is an exact p — 1 form for all formal solutions of the hydrodynamic system, since 
thereby /3 = dy is identically zero due to d'^S(t) = 0. 

The density p-form a and the flux p — 1-form (3 in the conserved integral fl2.19p satisfy a 
transport equation that arises from converting the boundary integral into a surface integral 




















( 2 . 20 ) 


through Stokes’ theorem and using the transport identity fl2.16p . This yields 

0 = — / cx + f d/3 = f T)tcy. + d/3 

dt Js{t) Js{t) Js{t) 

which holds on an arbitrary moving surface S{t) iff the integrand p-form vanishes, 

S)fa + d/3 = 0 . ( 2 . 21 ) 

Conversely, integration of this transport equation fl 2 . 2 ip over any moving surface S{t) yields 
a conserved integral (I2.19p . 

Note that if the conserved integral 02.191) is extended to the case p = n, with S{t) thereby 
being a moving domain V(f), then it coincides with the integral continuity equation 02.15P 
such that a = Te and /3 = <hje. 

2.3. Trivial conservation laws. A conserved integral (12.191) on a moving submanifold 
S{t), with any dimension 1 < p < n, reduces to a boundary integral iff the conserved density 
cx. = d0 is an exact p-form, holding for all formal solutions of the hydrodynamic system, 
where the p — 1 -form 0 is some function of the hydrodynamic variables and their spatial 
derivatives, and the time and space coordinates t, x. The corresponding flux is given by the 
p — 1-form /3 = —S)f0 from the transport equation fl2.2ip . If this flux /3 is non-zero then 
the resulting boundary integral has no physical significance, since the conservation equation 
for the integral is just an identity, 

^ [ d© = k / © = / s,© (2.22) 

Js{t) Jos{t) JdS(t) 

which follows from Stokes’ theorem combined with the integral transport identity (I2.16p . In 
this case the conserved integral fl2.19p and the corresponding local conservation law fl2.2ip 
are called trivial. 

However, if the flux in a conserved boundary integral fl2.22p is zero, 

'Dt& = Dt@ + Cu& = 0, (2.23) 

then the boundary integral is a constant of motion on the moving surface dS{t) of dimension 
p — 1, assuming dS{t) ^ 0. In this case the boundary integral itself is non-trivial, corre¬ 
sponding to a non-trivial local conservation law fl2.23p with a p — 1-form conserved density 
0 and with a vanishing p — 2 -form flux. 

When p = n, a trivial local conservation law on a moving domain is equivalent to a 
conserved density given by T = DivQ in terms of the vector function 0 = eJ0. The 
corresponding flux is given by $ = —DtQ — {Div Q)u. 

2.4. Determining equations. Necessary and sufficient equations will now be derived to 
determine all conserved integrals on moving domains and moving surfaces for the Euler 
equations fl2.4p - fl2.7p of inviscid compressible fluid flow. 

For fluid flow in an n-dimensional manifold M, a scalar function T will be a density for a 
conserved integral fl2.15p on a moving domain in the fluid iff DtT+{div u)T = D^T+Div {Tu) 
is a total covariant divergence —Div <h for some vector function $, where T and <1> depend on 
the time and space coordinates t, x, the fluid variables u, p. S', and their spatial derivatives. 
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Hence the defining equation for T and $ to be, respectively, a conserved density and a 
moving flux is simply 

DtT =-Div{Tu + ^). (2.24) 

The following result based on the variational bi-complex [15] gives necessary and sufficient 
conditions to determine T. 

Lemma 2.1. Let v be a tensor field on a Riemannian manifold M, and let denote 

the mth order covariant derivatives of v. A scalar function f{x, v, Vn,..., V^n) is a total 
covariant divergence Div F{x, v, Vv ,..., V^n) ijf 

EM) = 0 (2.25) 

where 

= + ( 2 . 26 ) 

m=l 

is the covariant spatial Euler operator (variational derivative) with respect to v. 

Here Grad^ denotes the m-fold product of the total gradient operator Grad; the super¬ 
script * denotes a formal adjoint defined by 

(7(6 ® ® e™, (GradMf) = i-MhiU ® ® 6 , Grad^f) (2.27) 

holding for arbitrary vector fields 6 ^md an arbitrary scalar function / on M. A proof of 
Lemma [ 2.11 employing index notation is given in Appendix B. 

Necessary and sufficient conditions for determining T are now obtained by applying this 
lemma to equation fl2.24p . 

Proposition 2 . 2 . All conserved densities T{t, x, u, p, S, Vm, Vp, VS", • • • 5 V^m, V^p, V^S) on 
a moving domain for the Euler eguations fl2.4p - fl2.6p of compressible fluid flow in an n- 
dimensional Riemannian manifold M are determined by the (necessary and sufficient) egua¬ 
tions 

K(AT) = 0, EflDtT) = 0, EsiDtT) = 0. (2.28) 

Moreover, a density will be non-trivial iff it satisfies at least one of the conditions 

EflT)fO, EflT) y^O, Es{T)fO. (2.29) 

This characterization of conserved densities has a straightforward extension to moving 
surfaces. 

A p-form function ck will be a density for a conserved integral fl2.19l) on a p-dimensional 
moving surface in the fluid iff 21 * 0 ; = DtCt -|- CuOt is an exact form —d/3 for some p — 1-form 
function /3, where o: and (3 depend on the time and space coordinates t, x, the fluid variables 
u, p, S, and their spatial derivatives. The following result based on the variational bi-complex 
[T5] gives necessary and sufficient conditions to determine o:. 

Lemma 2.3. Let v be a tensor field on a submanifold of a Riemannian manifold M 
with dimension n, and let denote the mth order covariant derivatives of v. A ho¬ 

mogeneous p-form function f{x,v,'Vv,...,'V^v) with l<p<n — lisan exact p-form 
dF{x, V, Vv ,..., V^n) iff 

d/ = 0 
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(2.30) 








where 


m=l 

is the total exterior derivative operator. 

Here A denotes the antisymmetric tensor product. A proof of Lemma [2.31 in local coordi¬ 
nates can be found in Ref. [12]. We remark that the homogeneity condition on / is necessary, 
as otherwise the p-form cohomology of M must be taken into account. 

From Lemma 12.31 a necessary and sufficient condition for determining ct is given by 
applying d to the transport equation 

Dtcy. + Cucx = -d/3 (2.32) 

which yields 

A(da) + Cuida) = 0 (2.33) 

since Lie derivatives and time derivatives commute with exterior derivatives. Then the Lie 
derivative identity 

Cuif) = u\df + d{u\f) (2.34) 

leads to a simple characterization of conserved p-form densities. 

Proposition 2.4. All conserved homogeneous p-form densities a.(t,x,u, p, S,'Vu,'Vp,'VS, 

..., V^u, V^p, V^S), with 1 < p < n — 1, on a p-dimensional moving surface for the Euler 
eguations fl2.4p - fl2.6p of compressible fluid flow in an n-dimensional Riemannian manifold 
M are determined by the (necessary and sufficient) eguation 

Dflda.) + d(MjdQ;) = 0. (2.35) 

Moreover, a homogeneous density will be non-trivial iff it satisfies the condition 

da ^ 0. (2.36) 

3. Main results 

We begin by recalling the notion of symmetries for Riemannian manifolds. 

A Riemannian manifold (M, g) possesses an isometry if there exists on M a vector field C 
satisfying the Killing equation 

C^g = 0. (3.1) 

From the Lie derivative identity Cc^g = V©C, the Killing equation is equivalent to V©C = 0. 
Here © denotes the symmetric tensor product. Similarly, a Riemannian manifold {M,g) 
possesses a homothety if there exists on M a vector held ( satisfying the homothetic Killing 
equation 

^(;g = = const. 7 ^ 0. (3.2) 

This equation is equivalent to V © C = Ap. 

A vector held y on M is curl-free (irrotational) if V A y = 0. Locally on M, this condition 
is equivalent to y = Vf, for some scalar held f. The identity 2VC = V©C + VAC shows 
that a curl-free vector held y is a Killing vector ( if and only if it is covariantly-constant, 
Vy = 0. 

We now state the main classihcation results for kinematic conserved densities on moving 
domains and moving surfaces in compressible huid how in Riemannian manifolds. The results 








are obtained by directly solving the respective determining equations in Propositions 12.21 and 
12.41 as carried out in section HI 


3.1. Conservation laws on moving domains. 


Theorem 3.1. (i) For compressible fluid flow fl2.4p - fl2.6p in a Riemannian manifold {M,g) 


of any dimension n> 1, the non-trivial kinematic conserved densities T(t, x, u, p, S) admitted 
for a general equation of state P{p, S) comprise a linear combination of 

mass T = p (3.3) 

volumetric entropy T = pf{S) (3.4) 

energy T = p{^g{u,u) + e) (3.5) 

(linear/angular) momentum T = pg{u,(), = 0 (3.6) 

Galilean momentum T = p{flj — tufS/fj), F-vpg = 0 (3.7) 

where e is the thermodynamic energy fl2.8p of the fluid, and f{S) is an arbitrary non-constant 
function, (ii) The only special equations of state P{p, S) for which extra kinematic conserved 
densities T{t,x,u, p, S) arise are the polytropic case 

P = a(5)pi+2/n 

with dimension-dependent exponent 7 = 1 + ^ where cr{S) is an arbitrary function, and the 
isobaric-entropy case 

P = k{S) (3.9) 

where k{S) is an arbitrary non-constant function. The extra admitted conserved densities 
consist of a linear combination of 

similarity energy T = p{g{u,f) — ]^Xt{g{u,u) +nP)), C^g = \g, VA = 0 (3.10) 

Galilean energy T = p{9 — tulVO + jXt‘^{g{u,u)nP)), Cyeg = Xg, VA = 0 (3.11) 
in the polytropic case fl3.8p . and 

non-isentropic (linear/angular) momentum T = pg{u,C)f{S), C(^g = 0 (3-12) 

non-isentropic energy T = Xpg{u,u)f{S) — / f{S)P'dS (3.13) 


in the is obaric-entropy case (IMP. where f{S) is an arbitrary non-constant function. 


The kinematic conserved integrals (12.151) corresponding to these conservation laws (13.3p - 
(I3.13P on an arbitrary spatial domain V(t) C M transported by the fluid are written out in 
Appendix A. 

The classiflcation presented in Theorem 13.11 generalizes a recent classiflcation [21 [3] of kine¬ 
matic conservation laws for the compressible fluid equations fl2.ip - fl2.3l) in R"" with equations 
of state P(p, S) that have an essential dependence on the pressure, Pp 7 ^ 0. In particular, 
the conserved integrals arising from the conserved densities fl3.3l) - fl3.1ip provide a covariant 
generalization (cf. equations (IA.ip - (lA.7p i of the well-known conserved integrals [9l [10] for 
mass, volumetric entropy, energy, linear and angular momentum, Galilean momentum, simi¬ 
larity energy and Galilean energy in R". The conserved integrals arising from the additional 
conserved densities (13.121) and (13.131) which hold for isobaric-entropy equations of state are 
apparently new. These two conserved integrals (cf. equations flA.8p and flA.Op ) are admitted 





































for any Riemannian manifold {M,g), including the flat case M"" (but they do not appear in 
the recent classihcation in M"" because equations of state with Pp = 0 were not considered). 

As a corollary of Theorem 13.11 note that there are no special dimensions n > 1 in which 
extra kinematic conserved densities are admitted. 


3.2. Conservation laws on moving snrfaces. 


Theorem 3.2. (i) For compressible fluid flow fl2.4p - fl2.6p in a Riemannian manifold {M,g) 
of any dimension n> 1, no non-trivial kinematic conserved p-form densities a{t,x,u, p, S) 
are admitted for a general equation of state P{p,S). (ii) The only special equations of 
state for which a non-trivial kinematic conserved p-form density cx.(t,x,u, p, S) arises is the 
barotropic case 

P = P(p). (3.14) 

The admitted conserved p-form density consists of 

circulation cx. = u (p = 1) (3.15) 

where u is the fluid velocity 1-form defined by the dual of u with respect to g (namely, 
Cjn = g{C,u) for an arbitrary vector field (). 


The corresponding kinematic conserved integral fl2.19p on an arbitrary curve (1- 
dimensional surface) S{t) <Z M transported by the fluid is given by the circulation 


d 

dt 



u = —{g ^{u, n) + e 



(3.16) 


where 

e(p) = j P~^P{p)dp (3.17) 

is the thermodynamic energy of the fluid. 

Unlike for conserved densities, no classihcation of kinematic p-form conservation laws for 
compressible huid equations have previously appeared in the literature. As a corollary of 
Theorem 13.21 there are no special dimensions n > 1 in which extra kinematic conserved 
1 -form densities are admitted, and no conserved p-form densities for 2 < p < n — 1 are 
admitted. 


3.3. Constants of motion. Finally, we state a classihcation of kinematic constants of 
motion on moving domains and moving surfaces in compressible huid how in Riemannian 
manifolds by examining when the net huxes in the kinematic conserved integrals vanish for 
all solutions of the huid equations. 

A moving domain V{t) <Z M necessarily has a non-empty boundary dV{t). Hence the net 
hux in a conserved integral fl2.15p on V{t) will vanish in general ih the hux vector $ itself is 
identically zero. In contrast, a p-dimensional moving surface S{t) <Z M (with 1 < p < n — 1) 
either has a boundary dS{t) or is boundaryless. If dS{t) is non-empty, then the net hux in a 
conserved integral fl2.19p on S{t) will vanish in general ih the hux p— 1-form (3 is identically 
zero, whereas if dS{t) is empty, then the net hux will always vanish identically. 

From the hux expressions in the kinematic conserved integrals given by Theorem 13.11 (cf. 
equations flA.lD - flA.9ll on moving domains) and Theorem 13.21 (cf. equation fl3.16p on moving 
curves), we immediately obtain the following result. 
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Theorem 3.3. For compressible fluid flow fl2.4p - fl2.6p in a Riemannian manifold {M,g) of 
dimension n > 1, the only non-trivial kinematic constants of motion are a linear combination 
of mass pdV and volumetric entropy pf{S)dV on moving domains V(t) C M, for 
any equation of state P{p, S), and circulation u on closed moving curves S(t) C M, for 
barotropic equations of state P{p). 

3.4. Hamiltonian symmetries and Casimirs. The well-known Hamiltonian formulation 
for the inviscid compressible Euler equations in M"’ (cf. [T2l ITT] ) has a straightforward co¬ 
variant generalization to an arbitrary Riemannian manifold {M,g). In covariant form, the 
Hamiltonian fluid operator is given by 

( {p~^curlu)\ —grad p~^gradS\ 

-div 0 0 • (3-18) 

— {p~^ grad S)\ 0 0 J 

This operator determines a Poisson bracket 

. /6G/6u\ 

hsF/6u 6F/6p 5F/5S)P \ 5G/5p \ dV (3.19) 

\6G/6S) 

satisfying (modulo divergence terms) antisymmetry and the Jacobi identity [I5], for arbitrary 
functionals P = f FdV and G = J GdV where F and G are functions of t, x, u, p, S, and 
covariant derivatives of u,p,S. Here 6/6u, 6/6p, 6/6S denote variational derivatives, which 
respectively coincide with the spatial Euler operators Eu, Ep, Es when acting on functions 
that do not contain time derivatives of u, p, S. 

The covariant Eulerian fluid equations fl2.4p - fl2.6p in (M, g) are given by 


/m\ /5E/6u\ 

dt[p\=P\ 5E/5p (3.20) 

\sl \6E/6Sj 

in terms of the energy density fl3.5p of the fluid. More generally, the covariant Hamiltonian 
operator Fi gives rise to an explicit mapping 

/6T/6u\ /u\ /r)“\ 

- H ST/Sp = X U = Uo (3.21) 

\6T/6S) \SJ \i,p 

which produces inhnitesimal symmetries (in evolutionary form) X = du J- g^dp J- g^ds of 
the compressible Euler equations fl3.20p from conserved densities T, where the components 
of the symmetry generator are given by 

ffl = —{p~^ curlu)\{5T / 5u) + grad{6T/6p) — {p~^gradS){6T/6S) 

g^ = — div {5T/5u) (3.22) 

f = {p-^gradS)\{6T/6u) 
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satisfying the symmetry determining equations 

Dtif' + g{u^ Grad grad )u — p~‘^f]'^grad P + p~^ Grad {Ppfj^ + PsV^) = 0, 

Dtp^ + {div u)r]^ + pDiv fj^ = 0, Dtp^ + g{ff‘., grad S) + g{u, Gradp^) = 0 

(3.23) 

for all solutions of the compressible Euler equations fl2.4p - fl2.6p . 

A conserved density T that lies in the kernel of the Hamiltonian operator fl3.18p determines 
a conserved integral called a Casimir [15] of the Hamiltonian structure. Every Casimir 
corresponds to a trivial symmetry, X = 0. From Theorem 13.11 a simple calculation shows 
that the only Casimirs arising from kinematic conserved densities T{t,x,u, p, S) are linear 
combinations of the mass pdV and the volumetric entropy pf{S)dV. In addition, 
all of the remaining kinematic conserved densities can be seen to give rise to non-trivial 
symmetries with components of the form 

- (xJV)m, pp = pp-Tpt-x\^p, = V^ -rSt-x\^S (3.24) 

which are equivalent to inhnitesimal point transformations X = rdt + xj + V^\ 9u + + 

p^ds on {t,x,u, p, S) where 77 “, pP, p^ are functions of t,x,u, p, S, while r, y are functions 
only of t, X. 

In particular, the kinematic conserved densities for energy fl3.5p . (linear/angular) momen¬ 
tum fl3.6p . and Galilean momentum fl3.7l) . which exist for a general equation of state fl2.7p . 
respectively yield the point symmetries 

time-translation X = dt, (3.25) 

isometry X = -|- ^g{u, curl()\du, (3.26) 

Galilean boost X. = t{grad'ip)\dx + grad'ip\du. (3.27) 

The extra kinematic conserved densities consisting of similarity energy fl3.10p and Galilean 
energy fl3.1ip . which exist only for a polytropic equation of state fl3.8p . and non-isentropic 
energy fl3.13p and non-isentropic momentum fl3.12l) . which exist only for an isobaric-entropy 
equation of state fl3.9p . yield the respective point symmetries 

similarity scaling X = Xtdt + — ^Xnpdp + ^{g{u, curl^) — Xu)\du, (3.28) 

Galilean dilation X = \Xt^dt + t{grad6)\dx — \Xntpdp + {grad6 — ^tXu)\du, (3.29) 

and 

generalized isometry X = f{S)^\dx — pC^f{S)dp + ^g{u, curl^)\du, (3.30) 

generalized time-translation X = f{S)dt + pCuf{S)dp. (3.31) 

Note the symmetry fl3.26p describes a space-translation symmetry if the Killing vector 
is curl-free (irrotational) and non-vanishing at every point in M, or a rotation symmetry 
if the Killing vector is not curl-free and vanishes at a single point (center of rotation) in 
M around which its integral curves are closed. If the Killing vector does not have either 
of these properties, then the physical interpretation of the symmetry (I3.26p depends on the 
nature of the zeroes and integral curves of C in M. 
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4. Solution of the determining equations 


In index notation flB.9|) - flB.lljl . the Euler equations fl2.4p - fl2.7p for inviscid compressible 
fluid flow in an n-dimensional Riemannian manifold M are written as 


- p-^D^P, 

(4.1) 

Pt = -Vi(pM*), 

(4.2) 

St = -u^WiS, 

(4.3) 

P = Pip,S), 

(4,4) 


where D^P = PpV*p + P 5 V*S'. 


4.1. Moving domains. A general kinematic conserved density has the form 


T{t, x',u\p, S). 


(4.5) 


Its total time derivative is given by 

DtT = Tt - + p-'(PpVV + PsV'5)) - TpV,(pR') - Tsu^VtS (4.6) 

on the space of (formal) solutions of the fluid equations fl4.1l) - fl4.4p . From Proposition l2.21 all 
kinematic conserved densities fl4.5l) are determined by the necessary and sufficient equations 


E^^{DtT) = 0, Pp(AT) = 0, Ps(AT)=0. (4.7) 


Expressions for the covariant spatial Euler operators E^i , Ep and Es are shown in index no¬ 
tation in equations fIB.lSD and flB.17D . The proper setting for evaluating E^i{DtT), Ep{DtT), 
Es{DtT) is the hrst-order jet space S) of the dynamical variables, which is coordi- 

natized by {t, x\ u\ p, S, 'Vju\ Vjp, V^P). Note that the covariant derivatives V* and V* will 
act on T only with respect to its explicit dependence on the coordinate x*. 

A straightforward calculation yields 


Ep{DtT) =Ttp + u^ViTp p ^PpV*T^i - pTppViE^ + p ^PpT^ijV^u^ 

+ (pPpTu^s ~ pPsTukp + PsTuk)p 

EsiD^T) =Tts + EV,Ts + p-^EsV^T^r - {pPpTu^s - pPsT^^p + P^T^Op'^VV 
+ p-^PsT^r^V^u^ + {Ts - pTps)V,u\ 

EuiiPtT) =Ttui + pViPp + u^VjT^i + pTppViP - p-^PpT^ijW^p - p-^PsT^i^V^S 
+ {pTps - Ts)V.S + (T„, - pT^^p)V,u^ + {pT^^ - 


(4.8) 

(4.9) 
(4.10) 





















By splitting each of these equations with respect to V^S', we get the system of 


determining equations 

Ttp + u^WiTp + p ^PpV^T^i = 0, ( 4 - 11 ) 

pPpT^ks - pPsTukp + PsTuk = 0, ( 4 . 12 ) 

P TppPjk — PpT^j^k = 0, ( 4 . 13 ) 

Tts + u^V^Ts + = 0, ( 4 . 14 ) 

P~^PsTj^k + gjk{Ts - pTps) = 0, ( 4 . 15 ) 

+ pVfcTp + u^VjT^k = 0, ( 4 . 16 ) 

9jk{T^i - pTuip) + 9ij{pT^kp - T^k) = 0, ( 4 . 17 ) 


which are to be solved for T(f, x*, m*, p, S) and P(p, S). We are interested only in non-trivial 
solutions, such that T and P each have some homogeneous dependence on at least one of m®, 
p, S. Note P can be determined only up to an arbitrary additive constant. 

In these equations fl4.11l) - fl4.17l) . note that t,x\u\p, S are regarded as independent vari¬ 
ables, while pjk is a function of x® such that ViPjk = 0. Hereafter we assume 

n > 1, (4.18) 

^ 0 or Tp ^ 0 or Ts ^ 0, (4.19) 

Ppj^OoiPs^ 0. (4.20) 


To proceed, we contract equation fl4.17|) with p®l, which yields 

= 0 (4.21) 

due to the condition fl4.18l) . By integrating equation fl4.2ip with respect to p hrst and 
next, we obtain 

T = pA{t, x\ u\ S) + x\ p, ^). (4.22) 

Equation fl4.12p then becomes 

= O' (4-23) 

We will return to this equation later, since it gives a case splitting. 

Next, we hnd equations 04.131) and 04.151) simplify to give 

9jkBp -|- PpA^j^k = 0, PjkBs + Ps^uiuk ~ 0) (4.24) 

with 

B = B- pBp. (4.25) 

By taking d^i of equation 04.241) and using condition O4.20p . we get A^ij^k = 0 which gives 

A = ^gjkU^u^A{t, X®, S) + u^Cj{t, x®, S) + Co{t, x®, S). (4.26) 

After we substitute this expression back into equation 04.241) . we have 

Bp + PpA = 0, Bs + PsA = 0. (4.27) 

Then by integrating this pair of equations with respect to p and S', we obtain 

B = -PA + f PAs dS + Po(t, X®). 
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(4.28) 
























With the use of 


PpAs = 0 (4.29) 

which follows from splitting equation fl4.23p with respect to u\ we now substitute equation 
fl4.28p into equation fl4.25p and integrate, giving 


B = pe A + / PAs dS + Bo{t, x^) 


(4.30) 


where e{p, S) is the thermodynamic energy fl2.8p dehned in terms of P{p, S). 

Then combining expressions fl4.30p . fl4.26p . fl4.22p . we see that the solution of the de¬ 
termining equations fl4.12p . fl4.13p . fl4.15p . fl4T7p . up to the case splitting fl4.23p . is given 
by 

T = + e)A + pu^Cj + pCo + j PAs dS + Bq. (4.31) 

Note the term BQ{t,x^) in this expression is a trivial conserved density (namely, it does not 
satisfy condition fl4.19p i. and hence we will put 


Bn — 0. 


(4.32) 


Substituting T from equations fl4.3ip and fl4.32p into the remaining determining equations 
fl4.1ip . fl4.14p . fl4.16p . and using equation fl4.29p . we get the system of equations 


2V(jCj) -|- A^pij — 0, 

Vii = 0, 

Cit + VjCo -l- {Pp pep pep)'ViA = 0, 

pCits P pVjCos -l- p{e'ViA)s P (PVtAjs = 0, 

Cot P PpV^Ci -|- (e pep)At = 0, 

pCost P Ps'^^Ci P p{eAt)s P PAts = 0, 

for A{t, x\ S'), Cj{t, x*. S'), Co{t, x*. S'). 

First, from equation (14.341) . we have 

A = Ao{t,S). 

Next, we contract equation f|4.33p with which gives 

WC, = -fioi. 

Substituting these expressions (I4.4np and (14.391) into the remaining equations in the system, 
we hnd that equations fl4.35l) and (14.361) reduce to a single equation 

Cu P V^Co = 0, (4.41) 

while equations (I4.37p and (14.381) become 

Cot P {pe — ^P)pAot = 0, (4.42) 

Cost P {es — fp ^.^ 5 )^ 0 * + (e + p ^P)Aots = 0. (4.43) 

We apply Vj to equation (I4.42p . which gives VjCot = 0. Integrating this equation, we get 

Co = F{x\S)pG{t,S). (4.44) 
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(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 


(4.39) 


(4.40) 

























































(4.45) 


Now we integrate equation fl4.4ip with respect to t, which yields 

Cj = -tVjF + Hj{x\S). 

Next we apply to equation fl4.40l) . giving Aotu = 0. Hence we have 

Ao = IoiS) + thiS)+t^h{S). (4.46) 

Then by integrating equation fl4.42l) with respect to t, we get 

G = (fP - pe)p{th + t‘^ 12 ) + J{S). (4.47) 

We note dp of expression (I4.47p yields (|P — pe)pp{tli + FI 2 ) = 0. After this equation is 
separated with respect to t, it is equivalent to the equation 

(P - lpe)ppA,, = 0. (4.48) 

Also, we hnd equation (I4.43p then reduces to give 

(p-^Ps - Pps + lpeps)Aot = (Pp - lipep + p-^P))A,,s. (4.49) 

These equations (I4.49p and (14.481) give case splittings, which we will return to later. Last, 
equation 04.331) separates with respect to t, yielding 

2V{iHj) = -Qijh, ViVjP = pijh- (4.50) 

The expressions 04.47p . 04.46p . 04.45p . fl4.44p . 04.39p . fl4.3ip . together with equation O4.50p . 
comprise the general solution of the determining equations 04.1ip - 04.17p . up to the case 
splittings 04.49p . 04.481) . 04.23p . 

Finally, we consider the various case splittings. From equation 04.48P combined with 
expression 02 . 8 p . we directly have 

Ppp - lp~^Pp = 0 or Aot = 0. (4.51) 

Similarly, from equation 04.23P combined with expressions 04.26p . 04.391) . 04.451) . we also have 

Pp = 0 or ios = Cks = 0. (4.52) 

These equations produce three distinct case splittings, as determined by the p dependence of 
P. The remaining case splitting is given by equation 04.49p . which has a more complicated 
form 

((1 + ^)p ^Ps ~ Pps)Aot = PpAots (4.53) 

depending on the p and S dependence of P. 

Case 1: Pp = 0. 

From this condition, we have that P is given by the equation of state P = P{S) and hence 
e = —p~^P{S) is the thermodynamic energy, where Pg 7 ^ 0 due to condition O4.20|) . Then 
equations 04.511) and 04.521) yield no conditions on Aq and Ck, while equation 04.53P reduces 
to the condition 

Aot = 0. (4.54) 

From the expression 04.46P for Aq, we see that equation 04.54p gives 

h = l2 = 0. (4.55) 


= 0, VtVjP = 0. 
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Then equation 04.501) becomes 


(4.56) 





































































This yields 


where 


V,F = Hi = 5;>'M(S)C(a 


)i 


.X 


(4.57) 


'^{iC{a)j) = 0, ViV ji}{a) = 0. (4.58) 

Since there are no further conditions, the expression fl4.3ip for the conserved density T 
becomes (after an integration by parts in the integral term) 


T = p{J + F)+ p{H, - tV,F)u^ + \phg,kU^v!^ - / hPs dS 


(4.59) 


with arbitrary functions /o(*S'), ^(*5'), p,(a){S), h'(^a){.S), and with arbitrary Killing vectors 
(C|^)(a;*), in addition to potentials ijj(^a){x^) for arbitrary curl-free Killing vectors. From the 
transport equation fl2.24p . a straightforward calculation now yields the moving flux associated 
to T, 

j {Hj - tVjF)Ps dS + u^ j IqPs dS. (4.60) 

Case 2: Pp ^ 0 and ^Ppp — p~^Pp = 0. 

By solving these conditions, we find that P is given by the equation of state P = a{S)p'^ + 
cTo(S'), with 7 = 1 and hence e = 7 ^cr(S')p'^“^ — p“^cTo(S') is the thermodynamic energy. 
Then equations fl4.5ip . (14.521) . fl4.53p yield 

^os = 0, Cks = 0, cTosAot = 0. (4.61) 

Substituting the expressions (I4.46p and (I4.45p for Aq and Ck respectively into equation (I4.6ip . 
we get 



Ji = Ji = const., I 2 = h = const.. 

(4.62) 


VjF = VjF{F), Hj = Hj{P)i 

(4.63) 

along with the 

case splitting 



(Jo = const. or Ji = J 2 = 0 . 

(4.64) 

If (To = const 

. then from equation (]4.5U]) combined with expressions (]4.63p. we 

have 


H, = + C,{P) 

(4.65) 

where 




= -gijhi h = 

(4.66) 


= ( 77 / 2 , h = 

(4.67) 


V(iCj) = 0 , ViVA = 0. 

(4.68) 


Since there are no further conditions, the expression (14.311) for the conserved density T in 
this case is given by 

T =p{J + F) + p{Hj - t'\/jF)u^ + + lit + Io){\pgjkU^u^ + ~ ^o) (4.69) 

with arbitrary constants Uo, /q, arbitrary functions o'(S'), J{S), and with an arbitrary homo- 
thetic Killing vector an arbitrary Killing vector in addition to a potential 6{x^) 

for an arbitrary curl-free homothetic Killing vector, and a potential 'ip{x^) for an arbitrary 
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curl-free Killing vector, where /i, I 2 are scaling constants associated to the homotheties. We 
note the term —( 12 ^^ + lit + Io)cfo in T is a trivial conserved density (namely, it does not 
satisfy condition fl4.19p h and hence we will put cxo = 0. Then, from the transport equation 
fl2.24p . the moving flux associated to T is given by 

¥ = {Hj - tV,F)g^^P + + lit + Jo)n'P. (4.70) 

Instead, if cxo 7 ^ const., then we have Ji = J 2 = 0. Consequently, equation fl4.50p combined 
with equation fl4.63p reduces to equation fl4.56p . which yields 

V,F = V,^(xO, H, = Cj{F) (4.71) 

where (j, 'ip satisfy equation fl4.68l) . Since there are no further conditions, in this case the 
expression fl4.3ip for the conserved density T becomes 

T = p{J + F)+ p{H, - tV,F)u^ + Io{\pgjkU^u^ + - uq) (4.72) 

with an arbitrary constant Jq, arbitrary functions ct(S'), cro(S'), J(S'), and with an arbitrary 
Killing vector (('■^ (x*), in addition to a potential '^(x*) for an arbitrary curl-free Killing vector. 
From the transport equation fl2.24p . the moving flux associated to T is given by 

¥ = {Hj - tVjF)g^^P + loU^P. (4.73) 

Case 3: Pp ^ 0 and ^Ppp — p~^Pp 7 ^ 0. 

In this case, P is given by a general equation of state P = P{p, S) other than the form 
arising in Case 2 , and hence the thermodynamic energy e has the general form fl2.8p . We 
then find equations fl4.5ip . (14.521) . fl4.53p yield the conditions 

^ot = 0, ^05 = 0, Cks = 0. (4.74) 

Substituting the expressions fl4.46|) and fl4.45|) for Aq and Ck into these conditions, we get 

h = h = 0, Jo = const., (4-75) 

VjF = VjF{P), Hj = Hj{P). (4.76) 

Consequently, equation (14.761) together with equation (I4.50p reduces to equations (I4.56P and 
(I4.7ip . Since there are no further conditions, the expression (I4.3ip for the conserved density 
T is therefore given by 

T = p{J + F)+ p{Hj - tVjF)u^ + Iop{\gjkU^u^ + e) (4.77) 

with an arbitrary constant Jq, an arbitrary function J(S'), and with an arbitrary Killing vector 
((’■^(x*), in addition to a potential V’(^*) for sji arbitrary curl-free Killing vector. Similarly to 
the previous case, the moving flux associated to T is given by 

¥ = {Hj - tvjF)g^^P + lovPP. (4.78) 

4.2. Moving surfaces. A general kinematic p-form conserved density, with 1 < p < n, has 
the form 

ai^...i^{t,x\u\p,S) (4.79) 

where Q;q...jp = a\i^...ip^ is totally antisymmetric. From Proposition 12.41 combined with the 
identity (IB.lSp . all kinematic conserved densities (I4.79P are determined by the necessary and 
sufficient equation 

T Pk{^ A (.P A ‘P)Pj]^ 0- 
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(4.80) 














































The proper setting for evaluating this equation is the second-order jet space J‘^{u\p,S) 
of the dynamical variables, which is coordinatized by p, SyV jU^V jPyV jS, 

V(jVjVfcP, Vj VfcS'), where, from relations flB.12p . 

V[jVfc]-u* =VpVqp = VyVqS'= 0. (4.81) 

Note that the covariant derivatives Vj and V* will act on only with respect to its 

explicit dependence on the coordinate x\ 

It will be useful to work with the dual of equation fl4.80p as follows. Let 

j-ji-j, ^ (4.82) 

which is a skew tensor given by the dual of with 

q = n — p, l<q<n (4.83) 

for convenience. Now apply equation fl4.8Up . yielding 


This is the necessary and sufficient determining equation for conserved tensor densities of 
rank q. 

The total divergence of is given by 

^ ^ VfcU* + pVkp + (4.85) 


A straightforward calculation of the terms in equation fl4.84p then yields 


+ 2VQ-|T^'i-4^-i'^5V|fc)5 + 

+ + 2T^'i-^-i%5V(fcpV,)^), 


(4.86) 


(4.87) 


(4.88) 
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and 


+ p-'(PpVV + PsV^S)) 

+ p-\PpV^p + Ps^'S)) 

- + p-\PpV^p + PsV^S)) ^ 

- + p-'(PpV*p + PsV^S)) (4.91) 

- T^’^-''^-^"ppVfcpV.(pn*) - T^’i-^-i%sn*VfcpV.5, 

- + p-^(PpVV + P^V'^)) (4.92) 

- r^'i-^-i'=5pVfc^V,(pn*) - 


on the space of (formal) solntions of the flnid eqnations fl4.1l) - fl4.4p . 

We now snbstitnte expressions fl4.86l) - fl4.92l) into the determining eqnation fl4.84p . combine 
terms after use of the derivative identities fl4.8ip . and split the resulting expression with 
respect to VjVkP, VjVkS, VjU^Vku"^, Vju^VkP, VjU^VkS, VjpVkP, VjSVkS, 

VjpVkS, 'VjU^, VkP, '^kS. This yields the system of determining equations 


= 0 , 

p^j.n-j,.pk^,gj)i ^ p^pn-j,.pk^,gj)i ^ 

{q - l)(5z[^'iTl^l-^-il\m + 

+ = 0 , 

(T^--^-^(\*(p-'Pp)p + = 0 , 

= 0 , 

(g - + pT^^-^'>-^\p5i^ 

+ + T=^-^<^-^\ipg^’‘p-^Pp = 0, 

(g - l)6i^nT\P-^,-i]k^ + 

+ 2Tn-j,-i[j ^Si^] + = 0 , 

(T^'i-^-i\.Psp + T^'^-'’^-^^„iPp)p-'g''' 

+ (T^--^-^^„.(p-'P5)p + T^"-^-^%..p-'P5)p*" = 0, 

- pVkT^^-^‘‘-^^pW + 

- (g - = 0, 

u^VjPn-ik^ _ ^^T^-d^.gikp-lp^ p pn-j,-ik^^ ^ 0 , 
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(4.93) 

(4.94) 

(4.95) 

(4.96) 

(4.97) 

(4.98) 

(4.99) 

(4.100) 

(4.101) 

(4.102) 














(4.103) 


u 








(4.104) 


which are to be solved for x\u\ p, S) and P{p, S). We are interested only in non¬ 

trivial solutions, such that T-^i 'I? and P each have some homogeneous dependence on at least 
one of M*, p, S, where P can be determined only up to an arbitrary additive constant. In 
these equations (I4.93p - (l4.104p . note that t, x*, u\ p, S are regarded as independent variables, 
while Qik is a function of x® such that V,g,i. = 0. Hereafter we assume the conditions fl4.18l) . 
dm, and (m, as before. 

To proceed, we contract equation fl4.93p with 6/ and integrate with respect to p, which 
yields 

pi-j, = S). (4.105) 

Then equation fl4.98p reduces to give 

= 0. (4.106) 

Next we combine equation fl4.1UUp with equations fl4.94p and fl4.1U5p to get 

SuiPp + fP-P-^\ip-^Ps = 0. (4.107) 

By taking d^i of this equation and using ds of equation fH.lOOp , we obtain 

fjJPs - pPsf) = 0. (4.108) 


From equations fl4.106p and fl4.108p we get = 0 due to condition fl4.20p . This 

gives 

fh-U = AP-P{t,x\ S) + u^BkP'P{t,x\ S). (4.109) 

Next, we see equation fl4.99p splits with respect to m®, yielding 

(g - = 0, (4.110) 

{q - = 0. (4.111) 


By contracting each of these equations with 6k , and using condition fl4.83p . we get 
and which gives 

j^h-A-d ^ X®), = BiP'P-^^{t, X®). (4.112) 

Then equation fl4.107p yields 

PsBiP'P = 0. (4.113) 

We will return to this equation later, since it gives a case splitting. 

We now note equation fl4.94p becomes = Q due to condition 04.201) . Hence 

(4.114) 


is a skew tensor. We then hnd equations 04.961) and 04.971) are identities, while equation 
04.95P becomes (g — l)5\iPB = 0. By contracting this equation 
with (5fc®®®, we obtain 

(n - 1 - q)BiP'P-P = 0. 

This equation gives a case splitting, which we will return to later. 
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(4.115) 














































Then combining expressions fl4.112p . fl4.109p . fl4.105p . we find that the solution of the 
determining equations fl4.93p - fl4.100p . up to the case splittings fl4.115p and fl4.113p . is given 
by 

Th-U = (4.116) 

where satisfies equation fl4.114l) . The term in T is a trivial con¬ 

served density (namely, it does not satisfy condition fl4.19l) ). and hence we will put 

= 0. (4.117) 

Substituting expressions fl4.116|) and fl4.117p for T into the remaining determining equa¬ 
tions fl4.10ip - fl4.104l) . and splitting with respect to u\ we get the system of equations 

= 0, (4.118) 

= 0, (4.119) 

- (g - -7 = 0, (4.120) 

= 0, (4.121) 

= 0, (4.122) 

V(, = 0. (4.123) 

From condition fl4.20p . we have that equations fl4.118|) and fl4.119|) yield V= 0. 
Then equation fl4.120p simplifies to give V= 0. This equation, combined with 
equation 04.1211) . yields 

(4.124) 

where is a covariantly constant skew tensor. Hence, equation 04.1221) becomes an 

identity, while equation fl4.123p holds as a consequence of the symmetries OB.161) of the 
Riemann tensor. 

Finally, we consider the case splittings. Clearly, we want Bi^^"'A ^ 0, otherwise T will 
be trivial. Hence equations fl4.113p and fl4.115|) directly give P 5 = 0 and q = n — 1. In 
this case, we have that P is given by the barotropic equation of state P = P{p) and hence 
the thermodynamic energy e has the general barotropic form e(p) = J p~‘^P{p)dp. Since 
g = n — 1 , the covariantly constant skew tensor P^F -in-i jg multiple of the volume tensor 

pF-P = Poe^'i-^T Po = const.. (4.125) 

Thus, expression 04.1161) becomes 

rph-B-i = (4.126) 

which is the general solution of the determining equations fl4.93p - fl4.lU4l) modulo trivial 
terms, where Pq is an arbitrary constant. 

From equation fl4.82p . the dual p-form density corresponding to the skew tensor density 
jpji-jn-i jg given by 

ai = Bogiju\ p = n- q = l. (4.127) 

By a straightforward calculation we find that the associated p — 1-form flux in the transport 
equation 02.321) is given by 

/3 = -BQ{\giju"u> + e - p~^P) 
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(4.128) 


























































which is a scalar (since p — 1 = 0). 


5. Concluding remarks 

The classification results in section |3] apply generally to inviscid compressible fiuids when 
the fiuid fiow is either isentropic (in which case S is constant throughout the fiuid domain 
M) or non-isentropic (in which case S is constant only along fiuid streamlines in M). This is 
a consequence of the case-splitting method that is used in section 0] to solve the determining 
equations. More specifically, although the isentropic fiuid equations fl2.ip - fl2.2p could pos¬ 
sibly admit additional conservation laws that do not hold for non-isentropic fiuid fiow, the 
determining equations show that this possibility does not occur. As a consequence, it turns 
out that all kinematic conservation laws of isentropic fiuid fiow arise from the kinematic 
conservation laws of non-isentropic fiuid by simply restricting the entropy S to be constant 
in M. 

It is worth emphasizing that these classification results are complete for conservation laws 
of kinematic form on moving domains and moving surfaces. Kinematic conservation laws are 
a physically important but limited class. In a subsequent paper, it is planned to classify all 
conservation laws of vorticity form on moving domains and moving surfaces. This class of 
conservation laws turns out to be much larger than the class of kinematic conservation laws, 
as shown by the examples of new conserved vorticity integrals found in recent work [1]. 

A fully complete classification of fiuid conservation laws will require going beyond those 
two forms for conserved densities and spatial fluxes, in particular by allowing dependence 
on arbitrary high order derivatives of all the fiuid variables. This open problem will remain 
a hard challenge for future work. 


Appendix A. 

Let V(t) C M be an arbitrary spatial domain transported by a compressible fiuid, and let 
0 be the outward unit normal on the domain boundary dV{t). 

For a general non-isentropic equation of state 02.71) . the kinematic conserved densities 
03.3p - 03.7l) yield the conserved integrals 


- pdV = 0 (A.l) 

Jv{t) 

L [ pf(S)dV = 0 (A.2) 

Jv{t) 

4/ p{l 9 {u,u) + e)dV = - [ Pg{u,i))dA (A.3) 

Jvit) JdV(t) 

4 / P 9 {uX)dV = - [ Pg{C,P)dA (A.4) 

Jv{t) JdV{t) 

4 f pi^lj - tu\V7p)dV = [ tPV^/J\PdA (A.5) 

dt Jv(t) Jdv(t) 
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where C(^g = 0 and = 0- the polytropic equation of state fl3.8p . the extra kinematic 
conserved densities flS.lOp and fl3.11l) yield the conserved integrals 


d 

dt 

A 

dt 





^Xt{g{u, u) + nP))dV 



Pg{i — tXu, i))dA 



p{9 — tuj V6' + \Xt^{g{u, u) + nP))dV 



tp{ye\p 


\Xtg{u^ i'))dA 


(A.6) 

(A.7) 


where C^g = Xg and Cyeg = Xg with VA = 0. The conserved integrals yielded by the extra 
kinematic conserved densities (I3.12p and (I3.13p for the isobaric-entropy equation of state 
fl3.9p are given by 


d 

dt 

A 

dt 



pg{uX)f{S)dV 



h{s)g{C,i')dA 



{lpg{u,u)f{S) - h{s))dV 


/ h{s)g{u, u)dA 

JdV(t) 


(A.8) 

(A.9) 


where 

h{S) = j f{S)P'dS. (A.IO) 

Note the conserved integral flA.4p describes linear momentum if the Killing vector C, is 
curl-free (irrotational) and non-vanishing at every point in M, or angular momentum if the 
Killing vector ( is not curl-free and vanishes at a single point (center of rotation) in M 
around which its integral curves are closed. If the Killing vector ( does not have either of 
these properties, then the conserved integral flA.4p can be viewed as describing a generalized 
momentum whose physical interpretation depends on the nature of the zeroes and integral 
curves of C in M. 


Appendix B. 

To begin, a complete transcription between geometric notation and tensorial index nota¬ 
tion will be listed. 


Notation. Vector product operations: 

a\b < — > a^bi, a A c < — > a(D c <—> 20 *- 

where a, c are arbitrary vector helds and b is an arbitrary covector held, on M. 
Tensor product operations: 


A\B i —^ 


if rankA > ranki? 
if ranki? > ra nk 4 


n.Jl-"3p 
J1 JP 


p\q\ p\q\ 


where A, B, C are arbitrary tensor helds on M. 
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(B.l) 


(B.2) 

(B.3) 


























Geometrical structures and operators: 

9 <—t Qij, e i —^ e ^^ (B.4) 

Riem <— > Rijk, R < —t Rikj^g^^ = R (B.5) 

V i—> Vi, div i—> Vi, Div <—> Di (B.6) 

grad i —^ g^^Wj = V\ Grad i —^ g^Wj = D\ (B.7) 

Grad"^ i — {Grad"^)* i—> ■ ■ ■ D^^)* = (-1)"^ (B.8) 

Fluid velocity, pressure gradient, and their covariant derivatives: 

u <—)■ M* (B.9) 

ujV — yu^Vi, divui — yViU^, curlui —)■ 2V^*m-^^, Vui — yViu\ (B.IO) 

GradPi—> D^P (B.ll) 

Covariant derivative identities: 

[Vi, Vj]/ = 0, [Vi, Vj]a^ = -Ri/a\ [Vi, ^ j\hk = Rijkk (B.12) 

where / is an arbitrary scalar function, a is an arbitrary vector field and b is an arbitrary 
covector field, on M. 

Lie derivative identities: 

Ru^ii-’-iq — U DkAi^...i^ + (3'^fc[i2---i5Vij]"?/ (B.13) 

= (g + + qD[i-^\{u'^Ak\i^...i^]) (B-14) 

= Dk{u^Ai^...i^) + (g + l)A[fci2...ig Vij]M^ (B.15) 

where A is an arbitrary skew tensor field, on M. 

Symmetries of Riemann tensor: 

^[ijk] 0, Ri[jk] 9lm Rm\jk\ 9li (B.16) 


Euler operator and its properties. The covariant spatial Euler operator fl2.26p is given 
by 

a.-+ (B.17) 


dv 

in the case of a scalar field v, and 
F • ■ ^ I 

ji-]q an*! J’lT" 


m>l 


dVkq ■ ■ ■ Vfc^n 


m>l 




d 


ki 


9V|i, 


(B,18) 


in the case of a tensor field In all cases, the Euler operator is uniquely determined 

by the following variational identity 

. U) + Dive[fM' -‘-n-i,) (B.19) 

3l"'3q 

holding for an arbitrary tensor field and an arbitrary scalar function 

f{x, '^kP^ '^^ji-jq, • • •, Vfci ■ ■ ■ where a prime denotes the Frechet 

derivative with respect to There is an explicit expression for 0 in terms of 

and partial derivatives of /, which we will not need here. The identity fIB.lQD 
leads to a simple proof of Lemma 12.11 
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If / = Div F then = Div and hence the identity flB.lOp 

holds only if (/) = 0 and = F'{w^^-Fj^,„-^). This 

3l‘‘‘3q 

directly implies {Div F) = 0, since the tensor held is arbitrary. 

Conversely, if E^i^...ip . . (/) = 0 then the identity fIB.lQD becomes f'{w^^"'^^ji-jq) = 

3l"'3q 

Div A homotopy integral can now be used to obtain /. Let 

be a one-parameter homotopy such that "^^jvjq 

with 1 ^( 0 )*^ being a hxed tensor held. Hence ^/| _ = 

= Dive{f\ dxV(^x)"^-^^h-jq), which implies f = Div F with 

F = Q{f\ _ , dxV(^x)^^ '^^ji-jq)dX -l- 00, where ©o is any vector held satisfying DwQq = 

f\ . 
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